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Abstract 

In this talk, we discuss how to estimate the gravitational lensing effect of a local void 
on the CMB polarization by using the LTB model. 

1 Introduction 

Type la supernova (SNIa) observations imply an acceleration of the cosmic expansion if the universe 
is homogeneous and isotropic on scales larger than 200Mpc. If we abandon this assumption called the 
Cosmological Principle, then other explanations become possible. The most interesting model of such 
a nature is the local void model, which was first proposed by Kenji Tomita in 2000 [1]. This model 
assumes that we are around the center of a low density spherically symmetric void and the spacetime is 
well described by the Lemaitre-Tolman-Bondi (LTB) model. In this model, the cosmic expansion rate 
decreases outward at each constant time slice, which produces an apparent acceleration of the universe 
when observed along the past light cone. Although this model violates the Cosmological Principle and 
requires the accidental situation concerning our location in the universe, it does not require any dark 
energy or a modification of gravity theory. Further, as far as the redshift-luminosity distance relation 
obtained by the SNIa observations is concerned, this model can reproduce the observational results with 
any accuracy because it contains at least one arbitrary function of the radius. Actually, it has passed all 
observational tests so far. Therefore, it is of crucial importance to find observational tests that enable us 
to discriminate this void model from the FLRW-based models, in order to establish the necessity of dark 
energy or a modification of gravity. 

One possible such test is to observe effects of the inhomogeneity on CMB temperature and polarization. 
For example, gravitational lensing is expected to generate B-modc and the (EB) correlation for an off- 
center observer in the local void model. In this talk, we explain how to calculate such gravitational 
lensing effects on CMB in the LTB model. 

2 CMB Polarization 

2.1 How to Represent Polarizations 

First of all, we explain the standard method to represent the polarization of radiations. Let us consider 
an quasi-monochromatic plane electromagnetic wave propagating toward an observer, and take an or- 
thonormal xy-ha&is that is orthogonal to the wave propagation direction. Then, the electric field of the 
wave is represented as E = E x e x + E y e y , with E x = a x sin(wi — e x ) and E y = a y sin(cji — e y ). 



In this setup, one can define parameters that represent polarization as follows: / := (a 



Q := (a y ) — (a x ), U := (2a y a x cos(e y — e x )), and V := (2a y a x sin(e y — e x )). These are called the Stokes 
parameters. Physically, / represents intensity (temperature), Q and U represent linear polarization, 
and V represents circular polarization. We ignore V because circular polarization is never generated by 
Thomson scattering in the early universe. 

If the orthonormal basis is rotated in the wave plane, Q and U are linearly transformed. Since this is 
inconvenient, we will introduce new quantities below that are independent of the choice of the orthonormal 
basis. 
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2.2 Polarization Distribution Patterns 

Up to this point, we have considered a wave propagating only in one direction, whose polarization can be 
described by the Stokes parameters (/, Q, U). In real observations, this set of parameters is measured for 
photons of each direction, and the result is represented by three functions on the sky, /(0 o bs)j Q(^obs) 
and U(6 ba), where o b s represents the position on the sky. 

In the present article, for simplicity, let us work in the flat-sky approximation. This approximation is 
valid when we consider only a small part of the whole sky. Let us define a tensor from Q and U as 



2 V U(x) -Q(aj) 

where the subscripts of P run over x and y (x = 1, y = 2). This tensor is called the polarization tensor. 

The polarization tensor field can be used to define the two functions on the sky, E{x) and B(x). 
which are independent of the choice of the orthonormal basis, by \7 2 E(x) = d a dbP a b(x), \7 2 B(x) = 
e ac dbd c P a b(x). Because these represent the 'gradient' and 'curl' (or 'rotation') components of the linear 
polarization distribution, respectively, they are called the E-modes and the B-modes, respectively. 

The power spectra and correlation functions are defined in terms of their Fourier transformations 
E(£) and B(t) as 

{XiMXtf)) = (2n) 2 S(£ + l')Cf^ , (2) 

with Xi, X 2 <E {0, E, B} and Xi, X 2 G {T, E, B}, where represents the intensity (temperature) 
fluctuation around the sky average ('2.7K'). If physics and the ensemble for averaging are invariant under 
a parity inversion, it turns out that Cj B = Cf B = 0. 

In the real, spherical-sky case, a similar argument holds. For details, the reader is referred to Ref. [2]. 



3 Gravitational Lensing Effects 

Inhomogeneous gravitational fields produce two effects on photon propagation. The first is a bending of 
its trajectory, and the second is the change of the photon energy in addition to the standard rcdshift by 
cosmic expansion. The latter is the so-called Sachs- Wolfe effect, which we do not consider in this article. 
Intuitively speaking, as far as CMB measurements by a fixed observer are concerned, the former — called 
'shear field effect' — can be further divided into two parts: (i) the change of the photon direction in the sky 
and (ii) the displacement of the intersection sphere of the past light cone and the last scattering surface 
in the direction perpendicular to this sphere. In order to give a definite meaning to this distinction, we 
need to introduce some reference FLRW model to define 'unperturbed' photon trajectories and past light 
cones. However, this procedure introduces the gauge freedom corresponding to the mapping between the 
real universe and the reference model, and thus make that distinction obscure. In fact, for the FLRW 
model with small perturbations, the displacement of the last scattering sphere can be set to be zero by 
an appropriate gauge choice, and in this gauge, the shear field effect can be represented only in terms 
of (z), namely, the 'gravitational lensing effect'. In the local void model, it is not so certain whether the 
same argument holds when the non-linearity of inhomogeneities is large. In the present article, we simply 
assume that the shift of the last scattering point in the direction normal to the last scattering sphere can 
be set to zero by a gauge choice. 



3.1 General Formula 

Under this assumption, the gravitational lensing effect on the CMB anisotropy can be simply determined 
by the two-dimensional shift vector d on the sky representing the difference between the observed direction 
of a photon and its initial direction on the last scattering sphere (Fig. la). The same formula holds for 
the LTB model as that for the FLRW model [3], [4], [5], and [6]. 

In the flat-sky approximation, let us represent the CMB temperature and polarization that would be 
observed in the reference exact FLRW model without gravitational lensing as 0, Q and U {E and B). 
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8<p(t s ) - 5ip(t Q ) 



Fig. la 



Fig. lb 



Figure 1: (a) Photon propagation in the LTB model. The black circle (dashed line) represents the last 
scattering sphere for the observer at the center O, and the red curve (solid line) represents that for an 
off-center observer P. The blue curve (arriving at P) is the photon trajectory, (b) The photon trajectory 
7 in the (r, (p) coordinates. Each trajectory is contained in the unique two-plane passing through O, P 
and the last scattering point. The green line (dashed-dotted line) represents the reference radial null 
geodesic 70. 



Then, the corresponding quantities observed with gravitational lensing are given by 

e(0 obs ) = e(8 ohs + d) = J -0^e l£ < e °^e(£) , (3) 

Q(8 ohs ) = Q(9 obs + d)= J -0^e^ e °^{E(£) cos 2^ - B(£) sin2<M , (4) 

U(8 ohs ) = U(9 ohs + d) = j ^^e^ e ^ + ^{E(e)sm2^ + B(£)cos2M , (5) 

where <pe is the angle between the vector t and the ^-direction on the sky. From this, the Fourier 
transformation of the differences are 

6m = I j0f®(£') w (£'> L ) . ( 6 ) 

SE(l) = J 0^{E(£')co S 2fa, e -B(£')sm2fa, e }W(£',L) , (7) 

5B(£) = J 0p{E(£')sm2ct)e, £ + B(£')cos24> £U }W(t',L) , (8) 
where W(£,L) = —£ ■ d(L), L =£ — £', <f>in := <p£> — <pg, and d(L) is the Fourier transformation of 

rf(^obs)- 



3.2 Null Geodesies in the LTB Model and the Shift Vector d 

Thus, the investigation of the gravitational lensing effect of a local void on CMB is reduced to determine 
d as a function of the photon direction. For that, we have to solve the null geodesic equation in the LTB 
model, whose metric can be written ds 2 = —dt 2 + S 2 dr 2 + R 2 (d8 2 + sin 2 8dip 2 ). Here R is a function of 
t and r, and S is written in terms of R and the curvature function k{r) as S — R'/(l — fc(r)r 2 ) 1 / 2 . 

In terms of the photon 4-momcntum = dx^/dX with affmc parameter A, the geodesic equation 
can be written as dp^/dX = —T^ p p y p p . Because of the spherical symmetry, this set of equations can be 
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reduced to the coupled ODEs for U), fi and p± denned by p* = u, p r = fiuj/S, and p\ = uj 2 (1 — ,u 2 ), where 
p± := R{(p e ) 2 + (p v ) 2 sin 2 d} 1 / 2 . Without loss of generality, we can assume that the photon propagate 
on the 2-plane with 8 = 7r/2, and therefore p e = (Fig. lb). Then, the geodesic equations are reduced 
to the set of four ODEs for u>(t), r(t), tp(t), and fi(t). 

It is impossible to express the general solution to this set of equations in terms of known functions. 
However, if we restrict the consideration to null geodesies passing through a point P close to the symmetry 
center O, we can find explicit expression for the solution in terms of integrals of known functions. Such 
a geodesic 7 stays close to some radial null geodesic 70. Hence, we can solve the geodesic equation 
perturbatively with respect to the deviation of the two geodesies 7 and 70. One subtle point in this 
perturbative approach is that 5(i turns out to show bad behavior near the observer. This problem can 
be avoided by introducing the variables b & c defined as b := r \J 1 — p 2 and c := r\i instead of r and fj,. 

The final result reads 



Sc(t) 
5u>(t) 



5c(t )exp (J^ 1 



(9) 



Sb(t) 



■ cxp 



Sb(t ) exp / dh I ^- - ^ + I + i 




SR R S cS 



dip(t ) ± 6b(t ) / dt 



to 



\c\R 



exp / dt 2 



R' R S 1 

1SR ~ R + S + cS 



(10) 

(11) 
(12) 



In Fig. lb, we can take the null geodesic passing through O in the direction f^bs as the reference 
geodesic 70. Then, it is easy to see that the values of 8b and Sc at present t = to can be expressed as 



Sb(t a ) 



Mo 



Sc(t ) = Dp 



(13) 



where D is the distance of the observer P from O and fXo = — cos(6* b s )- Then, (12) determines the shift 



vector d = {dg 



d v ), dg 



0, d v = S(p(t s ) — 5(p(t ), where t s is the last scattering time. 



4 Summary and Future Work 

In this paper, we developed a formulation to calculate the gravitational lensing effects on the CMB 
temperature and polarization for an off-center observer in a spherically symmetric void described by the 
LTB model. Explicit estimations of these effects are under investigation [7]. 
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